Current reversal with type-I intermittency in deterministic inertia ratchets 
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The intermittency is investigated when the current reversal occurs in a deterministic inertia 
ratchet system. To determine which type the intermittency belongs to, we obtain the return map 
of velocities of particle using stroboscopic recording, and numerically calculate the distribution of 
average laminar length (I). The distribution follows the scaling law of (I) oc e -1 ^ 2 , the characteristic 
relation of type-I intermittency. 
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In the last ten years there has been an increas- 
ing interest on the ratchet system. The ratchet sys- 
tem has been studied theoretically and experimentally 
in many different fields of science, e.g., molecular mo- 
tors Q, @, new methods of separation of particles 0, 
and current-voltage rectification in asymmetric supercon- 
ducting quantum interference devices More recently 
there have been several works on ratchet system in quan- 
tum domain 0,0, Josephson junction array 0, and etc. 

Ratchet system [g| is generally defined as a system 
that is able to transport particles in a periodic structure 
with nonzero macroscopic velocity, although on average 
no macroscopic force is acting For directional mo- 
tion of particle with no macroscopic force, or unbiased 
fluctuation in a periodic structure, the system has to be 
driven away from thermal equilibrium by an additional 
deterministic [Tol ITlj or stochastic pertubation [l2L Il3| . 
Besides the breaking of thermal equilibrium, for direc- 
tional motion the breaking of sp atial inversion symmetry 
is usually required further |14j . For this, the potential 
has a spatially asymmetric shape that is called 'ratchet'. 
It has also shown that the directional motion can exist 
in the presence of spatially symmetric potential with ex- 
ternal pertubation that is time-asymmetric [la ]. 

Many works on ratchet system have been mainly lim- 
ited on the overdamped case. In these works, systems in 
interest are related to microscopic scale in which thermal 
fluctuations or noises play a dominant role. Recently, 
Jung, Kissner, and Hanggi have considered the effect of 
finite inertia on ratchet system |16| . By considering the 
inertia term, the dynamics can exhibit both regular and 
chaotic behaviors. They have shown that deterministic 
chaos, to some extent, mimics the role of noise. Also, it 
shows that there exists multiple current reversals as the 
amplitude of external driving is varied. Thereafter this 
system has been called 'deterministic inertia ratchets' or 
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'deterministic underdamped ratchets' [r?l Hsj. 

In a recent letter, Mateos has shown that the origin of 
current reversal in the deterministic inertia ratchets may 
be related to a bifurcation from chaotic to regular regime 
[T^ | . The observation of intermittent chaos has also been 
addressed before the occurrence of current reversal. How- 
ever, Mateos has merely mentioned the diffusion property 
in intermittent chaotic regime. After his work, it has 
been conjectured that the mechanism of current reversal 
may be related to a crisis bifurcation in which a chaotic 
state suddenly becomes periodic 0]. Furthermore, it 
was shown that current reversal in the same system can 
occur even in the absence of bifurcation from chaotic to 
regular regime on the other parameter ranges. This cur- 
rent reversal is interpreted in terms of change of the sta- 
bility of the periodic rotating orbits At present, up 
to our knowledge, any further work on the classification 
of the type of intermittency just past the bifurcation has 
not been reported. 

On the other hand, the Pameau and Manneville types 
of intermittency are mainly classfied into types I, II, and 
III by the structure of local Poincare map, v n +i = v n + 
av^ + e, v n+ i = (l + e)v n + ai%, and v n+1 = —(l+e)v n — 
av^, respectively |'20j. And these types of intermittency 
are characterized by characteristic relations, (I) oc e^ 1 / 2 
for type-I, and (I) a r 1 for type-II and type-Ill, where 
(I) is the average laminar length. Here the parameter 
e in type-I intermittency is the channel width between 
the diagonal line and the local Poincare map, while 1 + e 
in type-II and type-Ill is the slope of the local Poincare 
map around the tangent point. 

The aim of this report is to investigate the intermit- 
tency of current reversal phenomenon in the determinis- 
tic inertia ratchets. We take the same system as in Ref. 
|19|. and explicitly show that there exists type-I inter- 
mittency when current reversal occurs from choatic to 
regular regime. 

Now, let us consider a system in which a particle moves 
in ratchet potential, subjected to time-periodic driving 
and damping. The equation of motion is written as 



d x , dx dV . . 

^ + 6 ^ + ^= aC ° S(w<) ' 



(1) 
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where b is the friction coefficient, and uj and a are the 
frequency and the amplitude of the external driving, re- 
spectively. Here V(x) in Fig. ^is the ratchet potential, 
and is given by 



V(x) = C 



sin27r(x — x ) + 0.25 sin47r(a; — x ) 
4^5 • 



(2) 



where C and xq are introduced to show the potential has 
minimum at x = with V(0) — 0, and 5 = sin(27r|a;o|) + 
sin(47r|xo|). 
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FIG. 1: The periodic asymmetric potential V(x) used, with 
parameters C = — (sin27ra;o + 0.25 sin 4ttxo)/4tt 2 5, xq = 
-0.19, and S~l. 614324 as in Ref. []]|. 

Among the three dimensionless parameters in Eq.Q, 
we vary the parameter a, and fix b = 0.1 and lo = 0.67 
throughout this paper as in Ref. |l9|. And the equation 
of motion, Eq.Q is nonlinear, and the inertial term as- 
sociated with the second derivative of x allows the possi- 
bility of chaotic orbits. We solve this system numerically, 
using fourth-order Runge-Kutta algorithms. Because of 
sensitivity of the above chaotic system, we calculate all 
values and parameters including 5 with double precision. 

In this system, there exists a bifurcation from chaotic 
to regular regime when current reversal takes place 
through varying the parameter a |l9j| . As shown in Fig. 
13 the particle shows intermittent chaotic behavior in the 
bifurcation region. The particle moves almost regularly 
in a negative direction, and occasionally shows chaotic 
burst. 

Furthermore, to study the behavior of the system 
during the current reversal, we plot a bifurcation dia- 
gram of velocity of a particle. For obtaining a bifurca- 
tion diagram, we take stroboscopic recording of the first 
derivative of x at times, t = fcr, where k is the posi- 
tive integer, and r is the period of the external driving, 
t = 2tt/oj. In the bifurcation diagram, Fig. |31 veloci- 
ties of the particle are plotted on the paramater range 
of a e (0.074000000, 0.086000000), after long initial tran- 
sient data are dropped. The result is analogous to Fig. 
2(a) in Ref. 0. 
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FIG. 2: For b = 0.1 and ui = 0.67, we show the intermittent 
chaotic trajectory of the particle at a = 0.080910000 during 
current reversal. 



0.35 




-0.25 1 1 1 1 1 

0.074 0.077 0.08 0.083 0.086 

a 

FIG. 3: The bifurcation diagram with varying the parameter 
a, other parameters are fixed, b = 0.1 and uj = 0.67. 

In Fig. |SJ two fixed points of velocities plotted in the 
bifurcation diagram at a = 0.074000000 correspond to 
the regular positive current of particle having period two, 
while four fixed points of velocities at a — 0.080990000 
correspond to the regular negative current of particle hav- 
ing period four\\sfy . During this current reversal, there is 
a period-doubling route to chaos, as shown in Fig. |3| Bi- 
furcation from chaotic to regular regimes takes place at 
critical value a c , just above a = 0.080947429 [2J. After 
this bifurcation, a periodic window corresponding to the 
regular negative current emerges. 

To determine which type above intermittency belongs 
to, we numerically obtain the return map, f 4 (v n ), which 
shows the relation between the velocity of particle, v n 
and the velocity of particle, v n+ \ after 4t time interval 
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elapse, with external driving period r, r = 2n/u. We 
use the value of r with double precision because of the 
sensitivity of the chaotic system. In Fig. 0] we plot the 
return map at just below a c . Fig. [S]is enlargement of the 
part of Fig. |2]in the vicinity of v n = —0.031. The points 
on the digonal line in the return map correspond to the 
states of v n +i — v n , i-e, four periodic motion. 
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FIG. 4: The return map, which is plotted at the parameter 
a = 0.080947429 just below a c . In the vicinity of -0.110, 
—0.031, 0.071, and 0.270, the return map is nearly tangent to 
the diagonal line. Four open circles indicate nearly tangent 
regions. 
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FIG. 5: The enlargement of Fig. 0Jin vicinity of v = —0.031. 
The return map seems nearly tangent to the diagonal line, 
but it is not exactly tangent. The channel width between the 
return map and the diagonal line is in the order of 0.000001. 



As shown in Fig. |SJ the return map is nearly tangent 
to the diagonal line at the parameter a just below a c . 
When a particle remains in a narrow channel, the trajec- 
tory shows almost a regular periodic behavior. After the 
particle escapes from this narrow channel, the trajectory 
shows a chaotic burst, and then goes back to the channel. 
(This process keeps repeating.) As increasing a from the 
value of just below a c to a c , the channel width narrows 
more and more so that the particle spends most of the 
time in the channel. After all, the return map crosses the 
diagonal line, when a becomes larger than a c , and two 
crossing points corresponding to the stable and unstable 
fixed points are made. Among these points, the stable 
fixed point corresponds to the regular negative current of 
particle. Like this, the intermittency emerges, before the 
return map undergoes a tangent bifurcation. This anal- 
ysis well agrees with the result shown in Fig. which is 
similar to that of Mateosfl9l|. 
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FIG. 6: The phase portrait of two attractors of the ratchet 
equation for b = 0.1 and u> = 0.67: one is the chaotic attractor 
for a — 0.080947429, just below a c and the other is the period 
four attractor for a = 0.080990000, represented by the center 
of four open circles. 

In Fig. one of two attractors is the chaotic one for 
a = 0.080947429, just below a c , and the other is the pe- 
riod four attractor for a — 0.080990000, corresponding 
to the regular negative current. Note that these two at- 
tractors are obtained by confining the dynamics in one 
potential well that include x = 0, i.e, in the range of 
x E (-0.37473443,0.62526557). The latter consists of 
four points at the phase space. In the chaotic attrac- 
tor, a particle spends most of the time in the vicinity of 
four point attractors corresponding to the regular nega- 
tive current, and once in a while intermittently moves in 
a chaotic way. 

So far, we have investigated the type of intermittency 
qualitatively. For quantitative characterization of the 
type of intermittency, let us survey the scaling law for 
duration of the laminar state, which means the particle 
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FIG. 7: The distribution of the average laminar length with 
varying the parameter a. It follows the scaling law of (I) oc 
e _1//2 with e = a — a c . 



remains in the narrow channel, as shown in Fig. [SJ We 



calculate the average laminar length (I) by averaging the 
iterated numbers of the return map at the laminar state, 
as varying a from the range of below a c to a c . Numeri- 
cally, we take the laminar state is the situation that the 
difference between v n+ \ and v is less than 0.0003. The 
result has shown in Fig. [7| The distribution of the aver- 
age laminar length follows the scaling law of (I) oc e -1 / 2 
with e = a — a c . It agrees with that the distribution of 
the average laminar length of the type-I intremittency is 
of the form {I) oc (T 1 ! 2 if there are no external noises 
[20II22I l23l| . Therefore, the intermittency that exists be- 
fore the bifurcation taking place from the chaotic to the 
regular regime is the type-I intermittency. 

In conclusion, we have investigated the mechanism of 
current reversal in deterministic inertia ratchets. By nu- 
merically obtaining the return map of velocities, and the 
scaling law giving the characteristic relation of intermit- 
tency, we have explicitly shown that the type-I intermit- 
tency exists when the current reversal occurs from the 
chaotic to the regular regime. 
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